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Abstract 

In this paper, we find various analytic (1+3)D solutions to relativistic ideal hy- 
drodynamic equations based on embedding of known low-dimensional scaling 
solutions. We first study a class of fiows with 2D Hubble Embedding, for which 
a single ordinary differential equation for the remaining velocity field can be de- 
rived. Using this equation, all solutions with transverse 2D Hubble embedding 
and power law ansatz for the remaining longitudinal velocity field will be found. 
Going beyond the power law ansatz, we further find a few solutions with trans- 
verse 2D Hubble embedding and nontrivial longitudinal velocity field. Finally 
we investigate general scaling fiows with each component of the velocity fields 
scaling independently, for which we also find all possible solutions. 
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1. Introduction 



Relativistic hydrodynamics [l|, Q is widely used in physics. Its particular 
application to high energy collision dynamics can be dated back to Laudau and 
Khalatnikov [3J in 1950's. For a brief history and an introduction to this aspect, 
see e.g. The relativistic ideal hydrodynamics (RIHD) has played a significant 
role in the past decade for the development of theory and phenomenology of 
the quark-gluon plasma (QGP) 0] — a deconfincd high temperature phase of 
QCD matter. The QGP can be experimentally created and studied by the 
Relativistic Heavy fon Collisions (RHIC). The space-time evolution at RHIC 
is characterized by fast longitudinal expansion (longitudinal flow) and strong 
transverse expansion (radial and elliptic flow). In particular, the large elliptic 
flow V2 in non-central collisions fij can be well described by RIHD [7] at low-to- 
intermediate transverse momenta for almost all particle species and for various 
centralities, beam energies and colliding nuclei, implying very low viscosity of 
the underlying matter, see reviews in e.g. @, [l^. It is now widely believed that 
the QGP in 1 — 2T^ region (as covered by current RHIC energy) is a "nearly 
perfect fluid" 11 1 l2| (see however caveats as discussed in e.g. l3| l3l)- The 



success of RIHD description for RHIC space-time evolution has become the basis 
for the paradi gm shift that the QGP in 1 — 2Tc region is actually a strongly- 
coupled Hquid 15, l3|. Such strong couphng and good fluidity close to Tc may 
have its origin in the magnetic component of QCD plasma [13, [3 which is 
ultimately connected with the mechanism of QCD deconflnement transition. 

The success of RIHD has also inspired a lot of interest in the formal study of 
RIHD. Analytic solutions from RIHD have been very useful for RHIC phe- 
nomenology, for example the well-known Hwa-Bjorken boost invariant flow 
[Til . [20I . There have been many efforts and progresses in finding analytic so- 
lutions, with an emphasis on solutions with potential importance for describ- 
ing RHIC collisions. Study aiming at (l-l-l)D solutions that give an alter- 
native description of the longitudinal expansion and a more realistic (non- 
boost-invariant) multiplicity distribution over rapidity can be found in e.g. 



21|, 122, |23|, 12J|. General Solutions of the (l-Hl)D RIHD with arbitrary ini- 
tial conditions in the forward light-cone for specific equation of states are also 
found in [2]|. Beyond (l-l-l)D, ellipsoidally expanding flows have been studied 
2^1 2^1, and generalization of the spherically symmetric Hubble flows to ac- 



celerating, non-Hubble flows in arbitrary (l-t-d)D has been done in 2l|. There 
have also been study of adding radial flow to a longitudinal Bjorken profile in 
e.g. [ill J ^^'^ more recently attempts to deal with elliptic flow analytically 
A new method for solving relativistic ideal hydrodynamics in (1-|-3)D is devel- 
oped in [2^1, where by explicitly embedding in the longitudinal and transverse 
radial flows, the authors reduce the RIHD equations to a single equation for the 
transverse velocity field only and find analytically all possible solutions with 
their transverse velocity fields having power dependence on proper time and 
transverse radius. 

In this paper, we aim to find a naly tic solutions to RIHD in (1+3)D, by taking 
the approach of Liao and Koch in 29| . Basically one uses lower-dimensional scal- 
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ing solutions (Hubble flow in ID and 2D) to build solutions in highcr-diinonsion. 
We will first study a class of flows with 2D Hubble Embedding, for which a sin- 
gle ordinary differential equation for the remaining velocity field will be derived 
in Scction-3. Using this equation, all solutions with transverse 2D Hubble em- 
bedding and power law ansatz for the remaining velocity field will be found in 
Section-4. More interesting solutions with transverse 2D Hubble embedding and 
nontrivial longitudinal flow beyond power law ansatz will be investigated and 
found in Section-5. In Section-6 we will study general scaling flows with each 
component of the velocity fields scaling independently, for which we also find all 
possible solutions. Finally we summarize and discuss the results in Section-7. 

2. Generalities 

We start with the relativistic hydrodynamics equations: 

T'"";„=0 (1) 

For ideal hydrodynamics, the energy-momentum tensor is given by 

j,mn = (e + -pg"''' = w u'^vT - p £f™" (2) 

with e,p and w = e + p being the energy density, pressure and the enthalpy den- 
sity defined in the flowing matter's local rest frame (L.R.F) which by definition 
are Lorentz scalars. The flow fleld u"^{x) is constrained by • Um = 1. In the 
usual {t,x) coordinates one can express u"^{x) as 7(1, u) with 7 = — -iP 

and V = dx/dt. 

To close the set of equations for flve independent fields e, p, v, we also need 
Equation of State (E.o.S) 

p = i'{e + p) = vw (3) 

The above means e = implying a speed of sound Cg = = 
in order to assure < 1, we require Q < u <l/2. 

2.1. Curved coordinates for 2D Hubble Embedding 

For a 2D Hubble Embedding study, we will use a coordinate system of 
{T±,r]±,(j)±, z): i.e. the transverse proper time, the spatial transverse rapid- 
ity, the transverse angle, and the longitudinal z. They are related to the usual 
{t, x, y, z) in the following way (with r± = \/x^ + y^) : 



rj_ = yJf^-T^ , r/_L = -In 



1, t + r^ 
-In 

2 t-rA_ 



and inversely 



1 , X + y ■ i 
2i X — y- 1 



t = T±_ cosh ?7j_ , z = z , 

X = T±_ sinh J7_L cos , y = t±_ sinh r}±_ sin (j)±_ (5) 
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The velocity field in these coordinates is related to u'^ = 'y{l,Vx,Vy,Vz) 
in flat coordinates {t,x,y^z) via 

u^-^ — 7[cosh77x — (vx cos4>_L + Vy sin0^) sinhjy^] (6) 
vP^ = — [(wa;COS0_L + Wj, sin (/)_L ) cosh — sinhry^] (7) 

u'''^ — (t;„ cos0_L — Wa; sin0_L) (8) 

T±_ smh ri±_ 

= -IV, (9) 

The metric tensor associated with the (r^, ri±, (f>±, z) coordinates is 

gran = Diag{l, -t\,-t\ sinh^ ^a, -1) 

g™=I?zag(l,-4,-^-^7^,-l) (10) 

smli r;^ 

With the metric above, it can be explicitly checked that u™Um — gmnU™"^" = 1. 

For the covariant derivatives in the hydrodynamic equations Eq.([T]), wc will 
need the Afline connections T'^^ = g^'Tkmn = g'^\{gkra,n + gkn.,m~ gmn.k)- 
In our case the non-vanishing connections are: 

_ _ -rV-L — r"?_L _ 
^ - ^-1 smn 7/^ , i , - i , - 



= -sinhT^^coshr;^ , F^^-,^ = r^|„, ^ (11) 



cosh ryi 
sinh7]j 



3. 2D Hubble Embedding and Derived Velocity Equation 

Now to find solutions in (H-3)D space-time, we first do an explicit 2D Hubble 
embedding, which involves the following ansatz for the transverse part of the 
velocity field: 

Vx = — = tanh77^ cos0^ , Vy = — = tanh77^ sm(f)± (12) 

which in the (t^, 77^, z) coordinates translate nicely into: 

u"^ = — ^ , m''^ = 

cosh r\\^ 

u-^-L = , , it^ = -iVz (13) 
The above can be further re-written as: 

7.™ = 7(l,0,0,z;,) (14) 
by introducing v, = v, cosh 77^ , 7 = -^/l — w^. 
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Now inserting the velocity field Eg. p^ together with the metric Ea.(fTU|) and 
the E.o.S Eq.([3]) (which substitutes the p by vw) into the energy-momentum 
tensor in Eq.(I5]), we can explicitly obtain all the non- vanishing components to 
be: 



rj-iT 



7 



7 «2 



rp<f>ct> 



P 



rj_ sinh 

-9-9 

, in 



(15) 

(16) 
(17) 

(18) 
(19) 



From the above, we can explicitly write down the hydrodynamic equations 
below: 

+ [27r_i (7,^^ +w^7,^ ) + 72 (2 + Tj^t;^ ,0] p = (20) 
P,ni- - (21) 
= (22) 
7^w^T_ip,^^ + [^'^vIt±_ + z/t_l] p,z 

+ [27^,(7,^^ +7,^ Vz)t^ + f{TA^Vz,r^ +21)2(1 + Ti_Vz,z ))] P = 0(23) 

The two equations about rj± and (p^ are trivially solved by setting all fields in- 
dependent of these two variables, i.e. p(x™) = p(t_l, z) and Vzix^) — Vz{t^, z). 
Then we can convert the remaining two equations into the following: 



(lnp),ri = J^[t^,z\ 



(24) 



-T±_Vz{Vz.r^ + VzVz.z) + V [2(1 - Z)^) -|- 2Tx_VzVz.r^ + r_i(l + vl)Vz.z] 



(-l-t-i/)i/ri(l-t;|)2 

(lnp),z = G[t±,z] 

_ T±{vz,T^ + VzVz,z) + v [~'2.Vz{l - vf) - rj_(l -h vl)\ 



(25) 

z.r^ - '2'T±VzVz,z] 



{-1 + V)VT^{1-VIY 

The commutativity of partial derivative (Inp)^,-^^^^ = (lnp)_2,Ti leads to the 
following constraint which is a derived equation only involving the velocity field: 



4 



+v[2vl - Urlvlv^^r^v,,, - 2vl{2 - 2t^v,,, + tI(v,\^ - 2v,,r^,z + v,^,)) 
-2v^{^l + 2t^v^,, + ri(3u^2^^ + 2w^,^^,^ + 3v^%)) + T±v^,{2v,^r^ + 

T±iVz,T^.T^ + Vz.z.z)) - T±iVz,T^i'2 + -ir^Vz.z) + T±{Vz^r^,r^ + Vz.z.z))] 

= (26) 

Once the above derived velocity equation has been solved, the matter field 
can be obtained as: 



p = Po- Exp 



Jt , n J Zn 



with po being the value at arbitrary reference point t±q , zq. 



(27) 



4. Analytic Solutions v^rith Transverse 2D Hubble Embedding and 
Longitudinal Power Law Ansatz 

As an application of the derived velocity equation (j26p in the previous sec- 
tion, we now find all solutions with power law dependence on the variables t± 
and z. Within such ansatz of the velocity field 

Vz = Arfz'^' (28) 

Eq. (P^ can be expressed explicitly as: 

(ACj.-AC2i.)r2+^z-2+c 

+ {-AB + AB^ + 2Av - ABv - AB'^v) t^z^ 

+ {A^B + 3A^B^ -AA^v- QA^B^v) t]^ z^^ 

+ [2A^v + A^Bu - A^B^v) rl^ z^^ 

+ {AA^BC - AA^Cv - %A^BCv) r^l^+^^z^^+^C 

+ (-A3C + iA'^C^ - (SA'C^y) r2+3B^-2+3C 

+ {AA^BC + AA^Cv - ^A^BCv) t\^^'' z-^+^^ 

+ {A^C + A^C^ - A^Cv - A^C'^v) r2+5S^-2+5C ^ q (29) 

In the above equation, terms with various powers of t^,z (and only power 
terms) appear: to make all of them, cither mutually cancel (among terms with 
exactly the same r^, p powers) or vanish by respective coefficients, to eventually 
zero is quite nontrivial. A thorough sorting of the sequences of , p powers can 
exhaust all possibilities to satisfy this algebraic equation. After enumerating all 
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the possibilities, we find only four sets of solutions: 

(i) A^O, with + < ; 

(ii) A = 1,B = -1,C = 1, with + + < ; 

{Hi) A = 1,B = 1,C = -1, with + < and + + > ; 

(iw) A = -3,B = -1,C = 1, with = ^ and x^ + + 9^^ < . 

(30) 

The solution (i) is just the 2D Hvibble expansion without z-niotion, the solution 
(ii) is actually 3D Hubble expansion, while the remaining two are more nontriv- 
ial. The applicable kinematic domain for each solution is given above in terms of 
(t, X, y, z) coordinates. The solution (iii) composed of two causally disconnected 
pieces: z > ^i"^ ~ ~JJ'^ and z < — -^/i^ — x^ — y^. It can be considered as an- 
other type of explosiorlfl, with the longitudinal flow being "anti-Hubble" . The 
solution (iv) corresponds to the special EOS v — 1/2, which will be further 
discussed in Section-6. 

5. Analytic Solutions with Transverse 2D Hubble Embedding and 
Nontrivial Longitudinal Flow 

In this Section, we try to seek more interesting solutions with nontrivial lon- 
gitudinal flow, i.e. beyond the simple power law ansatz. Instead of directly going 
to the second order differential equation in Eg. ipS)) . we return to Eas. ([55)l (|66 p 
noting the fact that both are homogeneous equations. Motivated by this fact 
and by the form of 2D and 3D Hubble flows, we start with the following joint 
ansatz for both the pressure and the velocity field: 

p{r^,z) = T-^g{-) (31) 

v.{Tj_,z)^f{ — ) (32) 

The above ansatz particularly features nontrivial dependence on the scaling 
variable ^ = and extra nonscaling time structure in the pressure. It is 
interesting to note that 2D Hubble flow corresponds to the trivial solution / = 
at a = 2 while the 3D Hubble flow corresponds to / = ^ at a = 3, both of which 
are special cases of the above ansatz. 

By inserting the above ansatz to (|63p and ([55]) , we obtain two first order or- 
dinary differential equations involving g{£^) = g{:f^) and /(^) = f{^)- Solving 
for f'lX) and g'{C) h'om the equations (with ' meaning ^), we obtain: 



•^Of course, the concepts of explosion and implosion depend on the reference point. We 
will choose the origin as the reference point if not further specified 
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.... ^ Kl-/(e)^)(2-a + a/(0^- 2^/(0) 

(34) 

We find that the Eq. (f33|) involves /(^) only and thus decouples from g(^). 
So one may first focus on solving and with /(^) obtained one can easily 
find g{^) from integrating Eq. p4|) . Again one may check that the 2D Hubble 
flow with / = at a = 2 and the 3D Hubble flow with / = ^ at a = 3 both 
satisfy the above equations. We observe Eg. ([55]) and Eq. fM]) are invariant un- 
der simultaneous transformations: f — *■ — f , / — *■ — /, g —f g. This looks like a 
manifestation of the time inversion symmetry, i.e. t± — > — t_l, in ideal hydro- 
dynamics, where dissipation terms are absent in stress energy tensor. However, 
we note tl > by definition. The symmetry can be interpretted as parity: 
z — *■ — It leads to the conclusion that any solution to Eq. lfSS]) and Eq. ([34|) is 
accompanied by its parity inversion counterpart. We will identify such pairs in 
explicit solutions later. 

To solve ([55]) . we can also consider it as (inversely) determining ^ = C(/) as 
a function of variable /. We can then introduce a new function 

.(/).^ (35) 
and recast ((33|) into a differential equation for t{f ): 

- [(^2 - a)f +{k^ + 2- a)t +{a- l)f t"^ + (a - 3)t^] = (36) 

In the above we have replaced the E.o.S parameter v hy k with the relation 
1/ = 73q7[('*^ ^ 1 such that f < 1/2). In what follows we try to seek solutions 
t{f) with two different types of ansatz. 

5.1. Solutions from Linear Polynomial Ansatz 

Motivated by the fact that all the coefficients in Eq. (|36|) are polynomials of 
both t and /, we naturally start with polynomial ansatz for the function t(f). 
The simplest form is: 

t{f)^af + p (37) 
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Inserting the ansatz into Eg. ([55]) . we obtain 



Co + CJ + C2f + C'sf = => (38) 
= Co = -/?[2 - a + + {2- a)a + {a- S)/?^] 

= Ci = a - + (-2 + 2a - K,^)a + (a - 2)a^ + (1 - a)(3^ + (9 - aa)^/?^ 
= C2 = -/3a[-4 + 3a + (3a - 9)a] 
= C3 = -a(l + a)[a + (a - 3)a] 

The resulting four algebraic equations can be easily solved, with all solutions 
listed below: 

1. a — 0, (3 — leading to 

t^O = ^ ' "^^^^ a = (39) 



II. a = K , (3 = leading to 



'=3^ ^ ^ = 2;^^ ' " = TT^ (40) 

The first solution coincides with the "anti-Hubble" solution found in previous 
section. The second solution is a new class of solutions. We will return to a 
detailed discussion of it at the end of this section. There are also nonphysical 
solutions which we discard, for example the solution a = — l,/3 = 0to the set 
of algebraic equations leads to t = —f and thus f — ±1 which does not comply 
with the physical constraint |/| < 1. 

One in principle can use various more general ansatz, for example: (a) poly- 
nomials beyond the linear one we used; (b) polynomial types with non-integer 
powers. These would typically be more complicated but may lead to more solu- 
tions. One class we examined is the form t — ap + /3 with 7 being any positive 
real number, which leads to no more solutions other than the ones found above 
with 7=1. 

5.2. Solutions from Nonlinear Ansatz 

In this subsection we go beyond the linear polynomial ansatz and show one 
example of nonlinear ansatz, which lead to interesting new solutions. Motivated 
by the 1- P term in Eq. ipS)) . we are led to the following nonlinear ansatz: 

tif) = af + P+^l^ (41) 

Furthermore we consider only specific solutions with a = 3 which has the ad- 
vantage of "killing" the term in the Eg. ([55)1 . Inserting this ansatz into the 
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equation for t{f) we obtain: 
1 



[Do + Dif + D^f + D^f + Z?4/' + D^f] = : 



(42) 



l_/2 

= Do = (/? + A)(l -K^ + a + p) 

= £)i = 3 - + (4 - K^)a + - 2/3(/3 + A) + (4 - K^)p + + 2ap 

= D2= P{k^ - 1 - 6a - 3p) + 3A(2 - a) 

= D3 = K^ -3 + {k^ -7)a- Aa^ + 20^ + 3p - 4ap 

= I?4 = 5q;/? 

= = 3a(l + a) 

By solving the above set of algebraic equations we found all the nontrivial and 

physical solutions listed below: 

III. a = 0, /? = 2, p = -2, A = -2 leading to 



t = -■ 



2/ 



/ = 



1 



1-/ •' 
IV. a = 0, /? = -2, p = -2, A = 2 leading to 

2/ , 1 



, with a = 3 and k = 1 



t ■ 



f 



l+f - 2+C 
V. a = 0, /? = l,p = —1, A = — 1 leading to 

/ 



, with a = 3 and k = 1 



1-/ •' 2(e-l) ' 
VI. a = 0, /? = —1, p = — 1, A = 1 leading to 

/ 



/ = ^ , with a = 3 and k = 2 



1 + / •' 2(e+l) ' 
VII. a = 0,/3 = 0,p=-l,A = leading to 



/ = ^ ^ , with a = 3 and k = 2 



t = - 



f 



l_/2 




i-C(0-C(C)] 

1 + i±i^c(c) + 

1 + i^c(o + '-^m 



with a = 3 and k = Vg 
27 

andC(0= ' 

m = 




16 



-i + rf-n/i-(i-T^a 



16 



(43) 



(44) 



(45) 



(46) 



(47) 
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The case for general a 7^ 3 with the same ansatz (HT|) is also investigated, in 
which case a set of 10 algebraic equations result. Many possibilities arise, and we 
analyzed most of them and found no more solutions, but few of the possibilities 
lead to really complicated situation that is not analytically tractable. Though an 
exhaustive investigation of general a is not achieved, it seems most likely there 
is no solutions other than those listed above with a = 3. There are certainly 
many more possible ansatz which could be tested in future works. 

5. 3. Discussion of Nontrivial Solutions 

In this section we choose to discuss properties and applicable kinematic 
domains of solutions II,III,IV,V,VI,VII. 

We start with solution-II. It is subject to the physical constraint |/| < 1 and 
reality condition of /. The applicable kinematic domain of the solution-II given 
in Eq. (|40p is given below (where we already replace ^ by the original coordinate 
variables z/t±^ and also replace by (1 — v)/v): 



l-3i 

constant x t. 



domain l:/^fZ^±V^SEMTZ^ 



2(1 -^.) 



with 2^/iy{l - v) < — < 1 (48) 



domain 2: f = 'J^^ V^V^i ~ Ml - 



2(1 - u) 



z 



with — < -1 (49) 



T 



2(1 

z 



with - 1 < — < -2yjv{\ - v) (50) 



domain 4: f^f/l^-v/^'T^?^^^^ 



z 



with — > 1 (51) 

In principle, any combination out of the four cases is an allowed solution. 
Here we will treat them as four independent solutions as there is no overlap 
in their applicable domains, i.e. they are causally disconnected from each other. 
Note solutions in domain 1 and domain 3 are related by parity inversion. The 
same is true for solutions in domain 2 and domain 4. The transverse part of the 
flow in all cases is always exploding 2D Hubble by embedding. The longitudinal 
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flow in z direction, as indicated by the sign of /, is also exploding outward in all 
four cases, since we have / > for > and / < for < 0. The pressure 
drops as tx increases keeping ^ constant. The solutions in some sense are like 
shock waves: at each given time moment r^, the flow is in certain domain of 
the full space; while with changing t_l, the spatial shape of the flow proflle (in 
z) does not change but only has its overall size scale up with tj_. 

Next we turn to solutions-Ill and IV. We first note they arc related by 
parity inversion. Therefore, we will concentrate on solution-Ill only. With 
physical constraint |/| < 1 applied, solution-Ill looks like: 



V. = ./( — 



-fi 1+1// 

p = constant x t_|_ 



(1//-1)^ 



with = i (52) 
domain 1: f = - — ^-i — 

with — > 3 (53) 

1 

domain 2: / = 



2 - z/t^ 



with — < 1 (54) 

Similar to solution-II, solution-Ill also contains two causally disconnected 
pieces. If we shift the scaling variable: z/t±_ z/t± — 2, which amounts 
to choosing z/t± = 2 as the reference point. The shifted solution / = — ^ 
represents imploding fluid. 

The other pair of solutions related by parity inversion is V and VI. We will 
focus on solution-V. In terms of original coordinate variables, the solution with 
physical constraint applied is shown as follows: 
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V. = /(-) 

p = constant x r , ^^^^ - — — ^ 

^ (1//- 1)15/8(1 + l//)5/8 

with y = - (55) 
5 

domain 1: / ^ 



1 + - 
2(^/ri - 1) 

1 - yS - Az/t^ 
2{z/t^ - 1) 



with — < -1 (56) 
domain 2: f 



with — < 1 (57) 

The transverse flow is again exploding by embedding, while the longitudinal 
part differs from solution-II. Note / does not change sign, solution in domain 1 
corresponds to fluid moving in the negative z direction, while solution in domain 
2 corresponds to fluid moving in the positive z direction. The pressure drops 
as T±_ increases keeping ^ constant. They can be interpretted as shock waves 
with infinite extension. 

Finally, we discuss solution- VII. Under the physical constraint |/| < 1 and 
reality condition, we are left with: 

p = constant x r'^^^il/f - 1)-^^ 

with y =\- 
7 

domain 1: / = = ^ [l - C(z/rx) - C(^/r±)] 

with — > or — < (58) 

domain 2: / = /2 = ^ 

with 0< — <lor -1< — <0 (59) 

Let us first comment on the reality of the solution: If 1 — (1 — f|^^)^ > 0, i.e. 
the square root in the definition of C and C is real, we can specify the argument of 
the square bracket such that C, and C, are complex conjugate to each other, which 
guarantees the reality of /i and /2. If the square root becomes imaginary, only 
/i remains real. The parity inversion simply maps /i and /2 to themselves. It is 
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also worth noting that solution in domain 1 and domain 2 contain discontinuity 
at = 0. They can be combined to give two smooth solutions at = 0: 

This solution is similar to solution-V. It can also be interpretted as shock 
wave with infinite extension. (|60p corresponds to fiuid moving in the positive z 
direction, while (|61|) corresponds to fluid moving in the negative z direction. 



6. General Scaling Solutions 

In this paper and in previous approaches, quite a few solutions have been 
found with simple scaling forms for velocity field, it is thus of interest to inves- 
tigate all possible solutions with the following simple scaling form velocity field 

(in t, X, y, z coordinates): 

X y z . 

Vx = oix- , Vy = ay- , = a^- (62) 

The three constants ax^y^z can in principle be different corresponding to anisotropic 
scaling solutions. These forms include many of the known solutions (cyclic in 
x,y,z): 

1) tta; = 1, ay^z = — Hwa-Bjorken or ID Hubble flow; 

2) ax.^y = 1, = — 2D Hubble flow; 

3) C(x.,y,z = 1 — 3D Hubble flow; 

4) setting az = 1 leads to a ID Hubble embedding in Liao and Koch paper (2^, 
in which a scaling solution ax.y = — 1 with v — 1/2 has been found; 

5) setting ax,y — 1 leads to a 2D Hubble embedding in the present paper, in 
which another scaling solution Uz = —3 with v = 1/2 has been found. 
Beyond the above, it should be mentioned that generalization of spherically 



symmetric scaling flow in arbitrary (l-|-d)D has been done in [21 1. 

It would be of great interest to exhaust all solutions with the above simple 
scaling ansatz. To do that, we simply submit the flow field into the hydrody- 
namics equation in ([1]) (and using simple flat coordinates [t,x)), which leads to 
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the following equations for p(t, x, y, z): 

[h [t^ - vh)] p^t + [htxax] p,x + [htyay] p,y + [htza^] p^^ 

+ [t [ast" - Cxalx" - Cyay - C.a'^')] P = (63) 
[htxax] p,t + [h {alx'^ + i/h)] p^^ + [ha^ayXy] p^y + [haxa^xz] p^^ 

+ [a^x {k.J^ - Kalx' - Ca^y' - ^^«'^')] P = (64) 
[htyay] p^t + [huyaxyx] p^x + [h {aly"^ + lyh)] p^y + [haya^yz] p^^ 

+ [%y {f^yt^ - ^t^lx^ - >^yC'ly'^ ~ C"^^^)] P = (65) 
[htzaz] p,t + [ha^a^zx] p^x + [ha^ayzy] p,y + [h {a^z^ + vh)\ p,^ 

+ \azZ (yKzt^ — C^CCx^"^ ~ ^t'^yV'^ ~ '^zCt1z'^)'\ p = (66) 

In the above, we have introduced a few constants (as certain combinations of 

^x,y,z^ ■ 

as = ax + ay + a^ (67) 

Cx,y,z = 2 + as - 2ax,y,z (68) 

K,y,z ^1 + as - ax.y,z (69) 

l^x,y,z = -1 + as + Olx.y,z (70) 
Cx.y^z ^ '^ax.yjZ i {^y^z^x ^z^x^y^ ('^-^) 

and also a coordinate-dependent function 

/i(x^) = - axx'^ - ayy^ - UzZ^ (72) 

The above derived equations (|63I64[ I65I66P can be recast into the following 
concise form: 

M^"" d^{\nY>) = B^" (73) 

from which we can explicitly solve out 

a,(lnp)=X;iB'' (74) 

The full expressions are somewhat long and not included here. In analogy to 
Section. 3, we can apply the compatibility condition, i.e. the commutativity of 
partial derivatives, which give six constraint equations on ax^ ay and az- Among 
these constraints, only four are independent. In practice, we find it easier to 
solve the constraint equations by starting from {lnp)^x,y = (^"•p),j/,x and then 
deal others subsequently. After tedious calculations, it has been found that an 
exhaustive list of all possible solutions include the following (all solutions valid 
subject to cyclic x ^ y ^ z ^ x): 

I. Qfj: = 1, ay^z = — Hwa-Bjorken or ID Hubble flow; 

II. ax,y = 1, = ~ 2D Hubble flow; 

III. ax,y.z = 1 — 3D Hubble flow; 
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IV. a,j; — ay and 1 + + cty + ctz — with i> = 1/2; 
'V. ax — and 1 + a^; + + az = with ly — 1/2. 

The appHcable kinematic region for all of the above solutions is simply a^x^ + 
a^y^ + a^z"^ < t^, which comes from the constraint that the flow velocity shall 
be less than the speed of light. We notice that there are no more solution cor- 
responding to general E.o.S parameter ly except the well-known Hubble flows in 
various dimensions. While at v = 1/2 (which somehow is a saturating extreme), 
we obtain two classes of solutions, with the solution IV including the solutions 
found previously as special cases and the solution V as a completely new class. 
Note the condition 1 + ax + ay + az — implies the matter is imploding(Q;i < 0) 
in at least one of the directions. The Solution IV and V represent implosion, 
but could contain one or two directions that are exploding(ai > 0). 

7. Summary 

In summary, we have investigated various analytic (1-|-3)D solutions to rela- 
tivistic ideal hydrodynamic equations based on known low-dimensional scaling 
solutions. We first studied solutions with transverse 2D Hubble Embedding, 
and derived a single ODE Eq.f^H) for the remaining longitudinal velocity field. 
Using this equation we found all solutions with transverse 2D Hubble embedding 
and power law ansatz (j28p for the remaining longitudinal velocity field, as listed 
in Eqs. ([5n)l . To find solutions with transverse 2D Hubble embedding but non- 
trivial longitudinal flows beyond power law ansatz, is both very interesting and 
extremely hard, and nevertheless we managed to find a few such solutions corre- 
sponding to different physical conditions, as listed in Eqs. (|39l40l43l44i I45I46I47|) . 
The last question we studied is all possible scaling solutions with the general 
scaling form in Eq. ([62|) for the velocity field, to which the complete answer was 
also found and listed. These analytic solutions describe various examples of ex- 
plosion, implosion, and shock-wave like processes in relativistic hydrodynamics. 
Though the initial motivation came from successful numerical hydrodynamic 
description of heavy ion collisions, the solutions presented in this paper may 
not directly apply to those but provide useful general indications for similar 
explosions in heavy ion collisions and also certain astrophysical processes. Ad- 
mittedly, a lot more work would be needed to reach a better and closer analytic 
description of the hydrodynamic phenomena relevant for today's relativistic 
heavy ion collisions. 
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